Impact of Step Defects on Surface States of Topological Insulators 
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The eigenstates in the presence of a step defect (SD) along xory axis on the surface of topological insulators 
are exactly solved. It is shown that unlike the electronic states in conventional metals, the topological surface 
states across the SD can produce spin rotations. The magnitudes of the spin rotations depend on the height 
and direction of the SD. The oscillations of local density of states (LDOS) are characterized by a wave vector 
connecting two points on the hexagonal constant-energy contour at higher energies. The period of the oscillation 
caused by the SD along y axis is V3 (-^) times that induced by the SD along x axis at a lai^ger positive (negative) 
bias voltage. With increasing the bias voltage, the period of the oscillation, insensitive to the strength of the 
SD, becomes smaller. At lower energies near the Fermi surface, the two types of wave vectors coexist in the 
LDOS modulations. These results are consistent qualitatively with recent observations of scanning tunneling 
microscopy. 

PACS numbers: 73.20.-r, 72.10.-d, 72.25.-b 
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Recently topological surface states have attracted much at- 
tention in the condensed matter community due to their poten- 
tial applications in quantum computing or spintronics [1,2]. 
The novel electronic states, which preserve time-reversal sym- 
metry, are produced by spin-orbit interactions. Such Dirac- 
cone-like surface states have been observed in three dimen- 
sional bulk insulating materials, such as Bii_tSbi [3], BiaSbs 
[4], Sb2Te3 [5], BizTea [5,6], TlBiSez and TlBiTej [7], by 
angle-resolved photoemission spectroscopy (ARPES). The 
surface energy band structure was determined by employing 
k ■ p theory [8], where an unconventional hexagonal warping 
term plays a crucial role in explaining the ARPES observa- 
tions. 

Scanning tunneling microscopy (STM) experiments have 
probed the electronic waves in the presence of step defect 
(SD) on the surface of the topological insulators Bi2Te3 [9,10] 
and the antimony (Sb) [11,12]. The absence of backscatter- 
ing of the topological surface states makes the local density 
of states (LDOS) near the SD more extraordinary as com- 
pared to that on the surface of conventional metals [13,14]. 
In Ref. [10], Alpichshev et al. observed the oscillations of the 
LDOS near a SD, dispersing with a wave vector that may re- 
sult from a hexagonal warping term. With increasing the bias 
voltage, the period of the LDOS modulation decreases. In this 
work, we investigate electron transport in the presence of a 
SD along x or y axis on the surface of topological insulators 
in the framework of quantum mechanics in order to explain 
the STM experiments. We treat the SD as a S{y) or S{x) po- 
tential barrier, similar to that in conventional metals [13,14]. 
We note that in Ref. [15], the authors studied the scattering 
from a 6{x) in strong topological insulators. However, they 
didn't take the hexagonal warping term into account, which is 
the key to understand the STM observations [9-12]. 

The momentum space Hamiltonian describing the surface 
states of topological insulators reads [8] 



= ( - + v{k,o-, - kyo-,) + A(f>{kj„ ky)o-,, (1) 



where / and cr,(; = x, y, z) are the 2x2 unit matrix and the Pauli 
matrices, respectively, m* is the effective mass of electrons, 
which is usually very large for the topological insulators, yu 
is the chemical potential, v is the strength of the Rashba spin- 
orbit coupling, the last term is the so called hexagonal warping 
term, and (piky^, ky) - k^ikl - 3^^). We note that the real space 
Hamiltonian corresponding to Eq. (1) can be obtained by tak- 
ing the transformations: kj^ = -/^ and A:,, - Obviously, 
the Hamiltonian (1) has the eigenenergies 



(2) 



with = and 1, and £k = yj A-cj)\k„ky) + v'^k'^ . 

We can see easily that x and y directions in the Hamilto- 
nian (1) are inequivalent due to the existence of the warping 
term. Therefore, it is naturally expected that different orienta- 
tion of SD leads to different oscillatory features in the LDOS. 
In the following we study two special SDs, which are usually 
observed in STM experiments. 

SD along y axis. We first consider the SD along y axis, 
which can be described by the S potential U(x) - Uq6{x) 
[13,14]. Here Uq is the strength of the SD and is usually weak, 
but rtfUo has a finite value. 

The wave function at the left of the SD, i.e. x < 0, is 



^s{ky, ky) 

rs{k„ky)\v(ik,-ky:) 



J{-k,x+kyy) 



U-kx,ky) 
rA-k„ky)\v(-ik,-ky) 



(3) 



w here ^k^k,) = A(P{k,,ky) + {-ly&k, r,(A:.„^,.) = 
^j^liky, ky) + v^k~, the first and second terms represent the in- 
coming and reflection wave functions, respectively, while the 
outcoming wave function at the right of the SD, i.e. x > 0, is 
found to have general form 
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FIG. 1: (Color online) The LDOS pi.ii(x, a>) as a function of distance 
from tfie line defect along y axis with different values of m'Uo and 
bias voltages. 



where C,, and are arbitrary constants to be determined by 
the boundary conditions at the SD, i.e. jc = 0, and 



{vkW-UKcl){k„K) + (-l)'A:,£k]. 



(5) 



We note that the first term in Eq. (4) is the tunneling wave 
function while the second term describes an extra spin rota- 
tion, distinguishing from the electronic wave function in con- 
ventional metals. It is nothing but this term that leads to the 
oscillations of the LDOS at x > [1 1]. 

Integrating the coupled Schrodinger equations associated 
with the Hamiltonian (1) in real space plus the 6{x) potential, 
we have the boundary conditions at the SD, 



Ys{k,,L)\v{ik,-ky) 



+e 



ik,Y I cos{k^x) + Bs sin(^,j-*^) 
Cs cos{k^x) + Ds sm{kxx) 



dx 



dx 



= 2m*f/ot/'/(0,3';k). (6) 



Substituting Eqs. (3) and (4) into Eq. (6), we obtain the spin 
(4) rotation constants 



J 



Ts(k,,k,)kAikx-m*Uo) 



kxl^sikx, ky) + Usv(ik, - ky)] + v^iik^ - 2m*Uo)v(ik, - ky) 



Ts{kx, ky)kx{ikx - m'Uo) 



VskMikx - ky) - (ik^ - 2m*Uo)[^s(kx, ky) + UsViik^ - ky)] 

V^^^(l + Mj + V^) H- UsTfs + iVsTJs 

I 



(7) 



where = ^s(-kx, ky)v(ikx - ky) + ^s(kx, ky)v(ikx + k^), and the 
reflection and tunneling coefficients 



Ts{-k„k,) 

Rs = ; - [Asviik, - ky) - CsUkx, ky)]. 



= 1 - — — ^ [Asviik, + ky) + CsU-kx, ky)] . (8) 

Obviously, when Uq - 0, we have Cj - Ds = Rs - and 
Ts = 1. 

In order to compare with the STM experiments, we calcu- 
late the LDOS near the SD, which can be expressed as 



A-v>0,*:,„.s 



(9) 



Here we emphasize that the formula (9) only considers the 
contributions of the topological surface states with ky > 
so that we can know clearly the LDOS modulations induced 
by the incoming and reflection wave functions or the out- 
coming wave function. We note that the LDOS observed by 
STM experiments should also include the contributions com- 
ing from those surface states with kx < 0, i.e. pstm(\x\,(jj) = 
Pi{x,cl)) + pii{x,co). In our following calculations, we use 
the physical parameters of the topological insulator BiiTes: 
A = 250.0 eV■A^ v = 2.55 eV■A^ and n = 0.334 eV [8,10]. 



3 



(a) to=100 meV 



(b) m=-^20 meV 



°<, O.O-' 




-0.2 -0.1 0,0 0.1 0.2 

k,(A') 
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C.S cos(A:,.}') + £)s ■&m{kyy) 



(11) 



where y[j = Ws(^.,Cj - ^y£)j), S., = w^ikyCs + ^.v^).?), vv.s = 
(ivk^y^^s{k^,ky). The spin rotation constants C., and the 
reflection coefficient K,, and the tunneUng coefficient 7~j are 
determined by the following constraints 



i/^;(ji:,0;k) = iA;^(x,0;k) 



FIG. 2: (Color online) The constant-energy contours of the surface 
state band at different energies. 



1 (dt//Ux,Q;k) d4r'{x,Q;k)] 



dy 



dy 



In Fig. 1, we present the LDOS with diff'erent values of m*UQ 
at high positive and negative energies. Obviously, the ampli- 
tude of the LDOS modulation near the SD depends strongly 
on the strength Uo of the 6 potential and the bias voltage cj. 
However, if oj is fixed, both period and phase of the oscil- 
lation keep unchanged with increasing Uq. When w = 100 
meV, the period Ty{cij) a; 58.0 Ain both sides of the SD. When 
oj = -120 meV, the period Ty(a)) ~ 42.0 A. At the positions 
with the same distance from the SD, the LDOS at these bias 
voltages has a maximum value and a minimum value, respec- 
tively. 

In order to understand the oscillatory features of the LDOS, 
we plot the constant-energy contours of the topological sur- 
face state band at w = 100 meV and -120 meV in Fig. 
2. We observe that the periods of the LDOS modulations at 
these energies are associated with a wave vector connecting 
two points on the corresponding constant-energy contours. In 
other words, r,(100) ~ ^ = ^ = ^ = 57.2 A, 

and r,(-120) ~ ^ = ^ = ^ . 43.2 A. There- 
fore, such oscillations of the LDOS are due to quasiparticle in- 
terference. It is obvious that there is no backscattering of the 
topological surface states in the LDOS, which is associated 
with the wave vector connecting two crossing points between 
the constant-energy contour and the x axis. 

SD along x axis. Now we consider the SD along x axis, cor- 
responding to that observed in BiaTea by the STM experiment 
[10]. Similarly, the wave function at y < has the form 



i/r;(x,y;k) 



e 



iikj(X-kyy) 



T,(k,,ky)\v{ik, + ky) 



(10) 



where we have used ^s(kx, -ky) = ^s(kx, ky) and Fsik^, -ky) = 
Ts{kx, ky). The wave function at y > is 



= UQilf]{x,0;k). 
Solving Eq. (12), we have 



(12) 



U-^,{kx, ky)af - U^viik, - ky)a 



,12 



,11„22 . „12„21 



r,(^,,^,.) 

C U^vjikx - ky)a'^ + U-Ukx,ky)af 
rAkx,ky)' al'af + al^al' 



^Ak^ky) 

K = mMikx - ky) - CsUkx, ky)], 



T,^l + l^^lM. [JlXik, + ky) - CsUkx, ky^] , (13) 

where = -j-^, C = ^s(kx, ky)[v{ik, - ky) ± v{ik, + ky)], 
and 



fl^ = ws^-ikl - \kxV') + (iky - ^U-)[wA.C - ^fsikx. k,)]. 



af = w,ky[^-(k, + if/-) - Csiiky - \v-)]. 



,21 



hj^Cs + (iky - ^u^)(c + 2w,kyk'), 



af = ky[(; + wsv\^{2iky - U^)]. 
Correspondingly, the LDOS in this case is 



..2;,2, 



(14) 



ilj]j{x,y;k)^T, 



„i{k,x+k,y) 



Ukx,ky) 
rAk„ky)\v(ik,-ky) 



p,,n{y,uy)^ ^ |tA;,„(x,y;k)|25(d^ - £k.,). (15) 

k„k,.>Q,s 
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FIG. 3: (Color online) The LDOS pijiiy, co) as a function of distance 
from the line defect along x axis with different values of Co and bias 
voltages. 
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FIG. 4: (Color online) The LDOS as a function of distance from the 
line defect along y oi x axis at zero bias voltage. 



smaller. When a) decreases to the values near the Dirac point, 
the quasiparticle interference associated with the wave vector 
- Qh or - qc becomes very weaker and the periods of 
the oscillations become very larger Therefore, the LDOS is 
almost constant. These results are consistent with the STM 
observations [10,11]. 

Because the modulation wave vector at w = 100 meV is 
different from that at a) = -120 meV, it is expected that the 
two modulation wave vectors compete at small energies. Fig. 
4 shows the LDOS at zero bias voltage for the two kinds of 
SDs. Obviously, the LDOS modulations cannot be fitted by a 
wave vector connecting two points on the Fermi surface [10]. 

In summary, we have investigated the impact of the SD 
along X or y axis on the surface states of topological insula- 
tors. We discover for the first time that there are spin rotations 
when the topological surface states move through the 5 poten- 
tial barrier. The oscillations of the LDOS near the SDs are in- 
duced by quasiparticle interference. This agrees qualitatively 
with the STM experiments. The period and phase of the oscil- 
lations are independent of the strength of SDs at high positive 
or negative bias voltage. But the amplitudes of the oscillations 
are sensitive to the strength of SDs and the bias voltages. We 
also find that the oscillations of the LDOS at high energies 
induced by the SD along y or x axis are associated with the 
same points on the constant-energy contours. Therefore, their 
periods have special relations, i.e. T^{\a)\) — y[3Tx{\(A) and 
= -^7".t(-|w|), where u is large. We hope that such 
relations could be verified by future STM experiments. 

This work was supported by the Texas Center for Super- 
conductivity at the University of Houston and by the Robert 
A. Welch Foundation under the Grant no. E-1411. 



According to Eq. (15) and using the parameters in BiaTe^, we 
also calculate the LDOS near the SD along x axis at different 
bias voltages and strengths of the S potential, shown in Fig. 3. 
We can see that the amplitude of the LDOS modulation also 
changes with and d. When u is fixed, the period and the 
phase of the oscillations are also independent of Uq. However, 
on the points symmetrical about the SD, the LDOS at w = 100 
meV has the same oscillatory features with that induced by 
the SD along y axis. In contrast, when oj - -120 meV, the 
LDOS has two maximum or minimum values. We note that 
r,.(100) ~ 33.0 Aand r,-(-120) ~ 74.0 A. The oscillatory 
characteristics are also produced by quasiparticle interference 
between two points on the constant-energy contours in Fig. 2, 
similar to the previous case. We find that ri(lOO) ~ i^^^^^^i = 
-J2_ = = 33.0 Awhile r^-(-120) ~ = ^ 

IqA-qfll 0.0952 ^ x\ j Iqf-qd \<Ie-1h\ 

= 74.98 A. With increasing u, Iq^ - qd and Iq^ - qd 
become longer, and so the periods of the oscillations become 
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